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ABSTRACT 

We compare models of line-driven winds from accretion discs and single spherical 
stars. We look at the problem of scaling mass-loss rates and velocities of stellar and 
disc winds with model parameters. We find that stellar and disc winds driven by 
radiation, within the CAK framework, are very similar as far as mass-loss rates and 
velocities are concerned. Thus we can use analytic results for stellar winds to rescale, 
in a first order approximation, numerical results for disc winds. We also show how the 
CAK stellar solutions change when we take into account effects of very low luminosities 
or line-driving force. 
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1 INTRODUCTION 



Radiation pressure has been long recognised as effective at 
powering mass-loss from luminous stars. It has been shown 
that stars with luminosities as low as 0.1 per cent of their 
Eddington limit can produce a powerful high velocity wind, 
when the transmission of radiation pressure to a flow via 
spectral line opacity is included. Studies of stellar winds 
driven by radiation pressure due to lines began more than 
two decades ago (Lucy & Solomon 1970; Castor, Abbott 
& Klein 1975, hereafter CAK). Later theoretical work by 
Friend & Abbott (1986) and Pauldrach, Puis & Kudritzki 
(1986; hereafter PPK) modified the CAK method and re- 
produced well empirically estimated time-averaged mass- 
loss rates and terminal velocities for main-sequence and 
evolved OB stars. After these successes of the modified 
CAK method, increasingly more attention has been paid 
to the instabilities inherent in the line-driving mechanism 
(e.g., Owocki, Castor & Rybicki 1988, here after OCR; Puis, 
Owocki & Fullerton 1993) and multi-dimensional aspects 
of stellar winds (e.g., Bjorkman & Cassinelli 1993; Owocki, 
Cranmer & Blondin 1994). 

A strong radiation field can also be produced by accre- 
tion discs, which are believed to be important and common 
components in many astrophysical objects, for example, cat- 
aclysmic variables (CVs), high mass young stellar objects 
(YSOs) and AGNs. Many of these objects show evidence of 
high velocity winds coming from the disk itself (e.g., Drew 
1997; Mundt & Ray 1994; Weymann et al. 1991). Small won- 



der then, that radiation pressure has been proposed to power 
mass-loss also from accretion discs. 

Despite a deep and quantitative understanding of 
radiation-driven winds in stars it has not been straightfor- 
ward to demonstrate that radiation can produce strong and 
fast winds from discs. The basic difficulty is the intrinsically 
two-dimensional, axisymmetric nature of disc winds. First 
studies of disc winds were focused on finding semi-analytic 
time- independent solutions. Such an approach required some 
simplifications and assumptions which had an essential in- 
fluence on the final outcome (e.g., Vitello & Shlosman 1988; 
Murray et al 1995). Nevertheless these studies illustrated 
in more detail the basic differences between driving a wind 
from a stellar photosphere and a disc photosphere. For ex- 
ample, Vitello & Shlosman showed that an increase of the 
vertical gravity component with height from the disc mid- 
plane prevents a disc from producing a wind unless the ver- 
tical component of the radiation pressure also increases with 
height. 

Numerical treatments of the disc wind problem have 
proven to be more successful as they allow solving of 
the multidimensional dynamical equations from first prin- 
ciples. Relatively early work by Icke (1980; 1981) using two- 
dimensional, time-dependent numerical simulations showed 
how radiation pressure due to electrons can produce a disc 
wind. More recently Pereyra, Kallman & Blondin (1997, 
see also Pereyra 1997) presented numerical calculations of 
the dimensional structure of line-driven disc winds for CVs. 
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However their spatial resolution of wind near a disc was too 
coarse to capture the wind structure in that critical region. 

Proga, Stone& Drew (1998, hereafter PSD) developed 
another numerical treatment of disc winds. PSD have identi- 
fied the inner and near disc as the important spatial domain 
and have therefore used a non-uniform (up to 200 x 200) 
grid to ensure that the subsonic acceleration zone near the 
disc surface is well sampled. 

PSD explored the impact upon the mass-loss rate and 
outflow geometry caused by varying the system luminos- 
ity and the radiation field geometry. A striking outcome of 
PSD's study has been the finding that winds driven from, 
and illuminated solely by, an accretion disc yield complex, 
unsteady outflow. In this case, time-independent quantities 
can be determined only after averaging over several flow 
timescales. On the other hand, if winds are illuminated by 
radiation dominated by the central star, then the disc yields 
steady outflow. PSD also found that the mass-loss rate is a 
strong function of the total luminosity, while the outflow ge- 
ometry is determined by the geometry of the radiation field. 
In particular, for high system luminosities, the disc mass- 
loss rate scales with the effective Eddington luminosity in a 
way similar to stellar mass loss. As the system luminosity 
decreases below a critical value (about twice the Eddington 
limit), the mass- loss rate decreases quickly to zero. 

PSD's calculations have been motivated by and de- 
signed for the case of winds from CVs. Their findings are in 
qualitative agreement with the kinematics of CV outflows 
inferred from spectroscopic observations. In a second paper 
(Drew, Proga & Stone 1998; hereafter DPS) they considered 
both the star and the disc as a source of mass. For model 
parameters suitable for a typical high mass YSO, where the 
stellar luminosity is ~ 2 orders of magnitude higher than 
the intrinsic disc luminosity, their radiation-driven wind disc 
model might explain the extreme mass-loss signatures of 
these objects. Most recently Oudmaijer et al. (1998) fol- 
lowed DPS and calculated a model for a B star surrounded 
by an optically thick disc. The resulting wind configuration 
could match with the wind structure in B[e] stars. In par- 
ticular, Oudmaijer et al.'s results show that the kinematic 
properties of the wind might explain the observed spectral 
features in HD 87643. 

In this paper we look at the problem of scaling disc 
mass-loss rates and terminal velocities with model param- 
eters. We will use results presented by PSD and DPS and 
use their approach to calculate some new results. Addition- 
ally we will review results for stellar winds to compare them 
with numerical results for disc winds. For that we will cal- 
culate some new stellar wind models for a parameter space 
not explored before. These new stellar calculations are for 
relatively low luminosities, like those of white dwarfs and/or 
a very low increase of radiation pressure due to lines as in 
highly ionized^] or low metallicity winds. We would like to 

* Abbott & Friend 1989 considered a related problem of the ef- 
fects of ionizing shocks in a stellar wind on the terminal velocity 



Table 1. Ranges of the model parameters which we have ex- 
plored. 



Parameter 


Range 




C's 


1.53 X 1CT 3 - 


1.38 x 10" 2 


a 


0.4 - 0.8 




Mmax 


1380 - 4400 




r D 


- 1.21 x 10" 


2 


r* 


- 1.21 x 10" 


■2 



point out that we treat our comparison between models of 
stellar and disc winds quite formally as our goal is to identify 
what and how model parameters control the mass-loss rate 
and velocity of disc winds. The resulting scaling relations 
will allow a better understanding of line-driving in general 
and also give some tools to estimate disc wind properties for 
many astrophysical applications before calculating detailed 
multidimensional time-dependent numerical models. 

We describe our numerical calculations in Section 2, 
present and compare our results in Section 3, and conclude 
with a discussion in Section 4. The Appendix summarizes 
the analytic results for M and the velocity law in the CAK 
framework. 



2 NUMERICAL METHODS 

To calculate wind models we adopt the 2.5-dimensional hy- 
drodynamical numerical method described by DPS (see also 
PSD). These models incorporate a non-rotating star, and 
a geometrically-thin and optically-thick accretion disc. The 
star and disc can both be a source of radiation and mass. 
PSD's formalism allows the stellar radiation to be included 
both as a direct contributor to the radiation force and as 
an indirect component via disc irradiation and re-emission. 
Each point on the disc and star is assumed to emit isotrop- 
ically. The model takes into account stellar gravity, gas 
pressure effects, and rotational and radiation forces. The 
gas in the wind is taken to be isothermal. In the dimen- 
sional form of the equations of motion for the wind model 
studied by PSD there are six model parameters: the stel- 
lar Eddington number, T» = akcGm » ^ ne disc Eddington 
number, Yd = 4 ^°q^j , (T* = xTd using the PSD x pa- 
rameter; all other symbols have their conventional mean- 
ing), the normalised sound speed, c4, and three parame- 
ters describing the force multiplier, a, k and r] max (see Ap- 
pendix). For parameters suitable for CVs (see PSD's ta- 
ble 1) the unit time, r = ^rl/(GM,) = 2.88s, the unit 

and mass-loss rate. Their approach was to assume that the radi- 
ation force is abruptly cut off at an adjustable distance from the 
star. 
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velocity vq = \J GM* /r» = 3017 km s 1 , the unit mass- 
loss rate, M = 8ircr*/a e = 2.62 x 1CT 5 M yr -1 and 
c' s = c s /v = 4.6 x 1CT 3 . 

Our calculations for stellar winds are as described in 
DPS with some changes outlined below. Because of the as- 
sumed spherical symmetry for stellar winds we use spheri- 
cally symmetric initial and boundary conditions. We set all 
initial velocity components to zero, except the radial compo- 
nent, iv, for which we use the CAK velocity (see equations 
A12 and A13). The density profile is given by hydrostatic 
equilibrium in the subsonic region and the CAK density pro- 
file (i.e., using the continuity equation, the CAK velocity law 
and mass-loss rate) in the supersonic region. The wind base 
density, pw, was chosen to sample the subsonic part of the 
flow and to allow initial transients to disappear on a short 
timescale. Typically, we use p;j, = 1CP 10 g cm~ 3 . The radia- 
tion force is recalculated using the method described in PSD 
for the case without a disc. Here we would like to mention 
that in PSD, DPS and Oudmaijer et al. (1998) the radiation 
force due to lines from a star was calculated using equation 
(A5), not PSD's equation (C2). Thus all previous calcula- 
tions of those authors were with a finite disc correction. 

To confirm and extend the PSD results, we present some 
additional disc models with a = 0.4, 0.6 and 0.8 calculated 
for various disc luminosities and with the stellar radiation 
switched off. The disc wind models are calculated as de- 
scribed in PSD. 

In numerical calculations with a finite line-force cut-off 
we used two approaches to saturate the force multiplier. The 
first approach we considered is simply: 

{kt~ a for t > tmin 
(1) 
M ma x for t t m i n 



where t min is chosen so kt m a in = M max . 

Our second, more physical approach, is that M(t) satu- 
rates gradually as the wind material becomes more optically 
thin (see also PSD). Specifically, we follow OCR who intro- 
duced an exponential cutoff in the line distribution which 
ensures the force multiplier saturates smoothly: 



For stellar winds, we calculated a few numerical test 
models to compare with analytic results summarized in the 
Appendix. We found good agreement between the numerical 
and analytic results - typically better than a few per cent. 



3 RESULTS 

Table 1 lists the parameter ranges explored in our numeri- 
cal simulations. Generally, P* and Yd control all the other 
parameters because they determine the wind photoioniza- 
tion structure for a given chemical composition. We thus 
need to explore only T» and Yd provided we calculate self- 
consistently the wind photoionization structure and subse- 
quently the radiation force. Unfortunately, this is not feasi- 
ble due to the wind's multidimensional nature and its po- 
tential time variability. However using models calculated for 
stellar winds and some basic physical arguments we may 
consider coupling between a, rj max , and k without detailed 
photoionization calculations. Specifically, we are guided in 
this matter by the results showed by Abbott (1982) and 
Stevens & Kallman (1990) and references therein. 

Using results of previous radiative force calculations 
for single stars and arguments from atomic physics, Gay ley 
(1995) argues that the maximum force multiplier, M max is 
relatively constant and depends on metallicity only. Gay- 
ley found that M max ~ a few 10 3 . Stevens & Kallman 
(1990) calculated the radiative force for a stellar wind ion- 
ized by an external X-ray source in massive X-ray binary 
systems. Their predicted maximum force multipier is also a 
few 10 and is nearly constant for low ionization parame- 
ter. For a very high photoionization parameter, M max de- 
creases to zero as all ions in a wind lose all electrons and 
there is a reducing contribution to the radiation pressure 
from lines. Stevens & Kallman showed how the force multi- 
plier decreases with the photoionization parameter. We will 
return to this point in Section 4. 

Table 2 summarizes our numerical results for stellar 
winds while Table 3 summarizes PSD's results and ours for 
disc winds with r, = 0. 



3.1 Dependence of wind properties on Y and M, 



M{t) = kt 



(1 ~\~ Tmax^) 



(l-a) 



T (l"«) 
' max 



(2) 



where T max — tr\ max and T] max is a parameter related to the 
most optically thick lines. Equation (2) shows the following 
limiting behaviour: 



Previous models of stellar winds driven by radiation force 
have explored a wide range of model parameters suitable 
for hot luminous stars. In the Appendix we briefly describe 
the basics of the CAK method, with the finite disc correc- 
tion (hereafter CAK&FD), and with explicit introduction of 
Mma X (hereafter CAK,FD&M mal ). We also summarize the 
analytic results for M* and the velocity law. 



lim M(t) 



T m ax — ► OO 



lim M(t) 



kt 
M, 



771 (IX. C 



(3) 
(4) 



In Fig. 1 we show (a) model mass-loss rates and (b) 
the terminal/typical velocities as functions of the Edding- 
ton number. All models on Fig. 1 are for a — 0.6, k = 
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Figure 1. Model mass-loss rates (upper panel) and terminal/typical velocities (lower panel) as functions of the Eddington number. The 
solid lines plot results for stellar models in the CAK case (equations A13 and A14), while the thin dotted lines plot results for stellar 
models in the CAK, FDUM m ax case (equations A25 and A26). The thick dotted lines plot results for stellar models in the CAK&FD case 
(equations A18, A19, and A20) while the crosses correspond to numerical results in this case. The asterisks connected by the triple-dot 
dashed lines correspond to numerical results for stellar models in the C AK,FDkM max ,e case. The open triangles represent disc results 
for models in the CAK&zFD case while the open circles, connected by the dashed lines, represent disc results in the CAK, FDIkM max . e 
case. On upper panel the thick solid vertical line is for F = 1/(1 + M max fFD,c)- The alternative ordinate on the right hand side of the 
upper panel is the dimensionless wind mass-loss rate parameter, M' = M / Mq, and on the lower panel it is the dimcnsionlcss velocity 
parameter, v' r = v t /vq. 



0.2 and M max = 4400. We plot analytic results for the 
stellar mass-loss rate, M, , and the terminal velocity, Voo 
in four cases: CAK, CAK&cFD, CAK,FDhM max , and 
CAK,FDkM maXie . The last case, CAK, FDSzM max , e is 
where the finite disc factor and OCR's exponential cutoff 
in the line distribution are included. 

For high T», the results shown on Fig. 1(a) confirm the 
well known CAK relation M* oc fy Q (the thin solid and 
thick dotted lines). Additionally our results confirm that the 
inclusion of the finite disc correction reduces M» by a factor 
of ~ 2 « (1 + ct) 1/a , for a ~ 0.6 (Friend & Abbott 1986 and 
PPK). This figure also reveals that the CAK scaling does 
not apply for very small F, if we take into account the fact 



that the force multiplier can not increase to arbitrarily high 
value (the thin dotted line and the asterisks). Our analysis 
presented in Appendix shows that the CAK scaling does not 
apply for 

r* < 1/(1 + (1 - a)f FD , c M max ). (5) 
The mass loss vanishes at 

r„ = l/(l + M max f FD , c ). (6) 
Previous studies by other authors of line driven stel- 
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Table 2. Summary of numerical results for stellar winds with a = 0.6, k=0.2 and M n 



4400. 



r* 


M, 




v r (10r„) 




(M Q yr" 


l ) 


(km s- 1 ) 


C! AK FDkrM ease 
o s\i \ , jp l^ocivimax ccibe 








ISiv i n -4 

iJ.U^ A A-KJ 


9.0 x 10" 


13 


2770 


6.03 x 10~ 4 


4.1 x 10" 


12 


4550 


1.21 x 10~ 3 


1.4 x 10" 


■11 


7640 


3.84 x 10~ 3 


9.6 x 10" 


■11 


11130 


1.21 x 10~ 2 


6.5 x 10" 


■10 


10750 


CAK,FD&M ma x,e case 








3.84 x 10~ 4 


4.2 x 10" 


■14 


2630 


6.03 x 10~ 4 


6.3 x 10" 


■13 


3970 


1.21 x 10~ 3 


5.5 x 10" 


12 


5660 


3.84 x 10~ 3 


6.6 x 10" 


■11 


7940 


1.21 x 10~ 2 


5.6 x 10- 


■10 


9180 



lar winds were for luminous stars for which the pres- 
ence of M ma x was not important. However for r, S; 
0.01 in the CAK,FD&iM max case, M, is lower and the 
M* vs. r„ relation is steeper than in the CAK&zFD 
case. For r. > 0.01, M, in the CAK,FDhM m ax 
case converges with M* in the CAK&iFD case. On the 
other hand, M* in the CAK, FD&iM m ax,e case is steeper 
than in the CAK, FD8zM m ax case. Additionaly, M» in 
the CAK, FD8zM m ax,e case converges with M* in the 
CAK&zFD case for higher F, than with M» in the 
CAK, FDSzM ma x case. 



Comparing results shown in Fig. 1 , we find that 
for the CAK, FDhMmax case v x remains sensitive to 
Mmax for higher r, than M*. Quantitatively «oo in the 
CAK, FDSzMmax case starts to deviate from Vao in the 
CAKhFD case at 



r* = 1/(1 + 2(1 - q)/fd, c M„ 



while M* starts to deviate at 



Fig. 1(a) also shows numerical results for the disc 
mass- loss rate, M D , from PSD (with a new point for 
Yd = 6.03 x 10~ 4 , see Table 3), and for disc wind mod- 
els with infinite M ma x, the open circles and triangles, re- 
spectively. PSD results are analogous to stellar results in 
the CAK, FD&zMmax,e case, while disc results with infi- 
nite Mm a x are analogous to stellar results for the CAK, FD 
case. Comparison between these analogue models demon- 
strates that the stellar mass-loss rate and the disc mass-loss 
rate depend on luminosity in a similar way for all lumi- 
nosities. The most convincing case is for disc wind models 
with infinite M ma x where disc mass-loss rate predictions al- 
most match the results for the stellar C AKSkFD case. In 
the CAK, FD&zMmax,e case, lower Md approaches M» in 
the stellar CAK, FDSzMmax and CAK&zFD cases as the 
Eddington number increases. 

Fig. 1(b) illustrates the behavior of with the Edding- 
ton number. The CAK theory predicts that «oo stays almost 
constant with F* for r, < 1 (see equation A13). Friend & 
Abbott (1986) showed that the finite disc factor introduces 
a dependence of Vao on T* but it is a weak dependence. Our 
calculations confirm these results and also confirm that t>oo 
is higher for the CAKSzFD case than for the CAK case 
by a factor of ~ 2 (Friend & Abbott 1986; PPK). However 
if we take into account M ma x, Voo decreases strongly with 
decreasing T» (see also equation A22). 



r* = i/(i + (i - q)/fd, c m„ 



higher than the former by a factor of ~ 2. This is unsur- 
prising, given that the terminal velocity of stellar winds is 
a global quantity determined in the whole acceleration zone 
whereas M is a local quantity determined at the critical 
point and conserved in steady state models (see also Abbott 
& Friend 1989). 

The characteristic velocity of disk winds does not de- 
pend on the Eddington number for infinite M max . However 
for Mmax = 4400, the characteristic velocity is a strong func- 
tion of Td - Thus the disk wind velocity and the stellar wind 
velocity are sensitive to the Eddington number in similar 
ways. The main difference is that disc winds react to finite 
Mmax for higher Eddington number than stellar winds. This 
difference indicates that M(t c ) is higher for discs than for 
stars. 



So far we showed that a finite M ma x reduces M» and 
Voo significantly unless M(t c ) < M max . In the regime where 
M(t c ) < Mmax the actual value of M ma x does not matter 
and M t and Vao do not depend on M m ax nor r/ m ax- As CAK 
found the main stellar properties are determined at the wind 
critical point. For example, equation (A18) shows how the 
mass-loss rate depends on the force multiplier at the criti- 



© 0000 RAS, MNRAS 000, 000-000 



6 D. Proga 



Table 3. Summary of parameter survey for disc winds with k = 0.2 and M„ 



■ 4400. 



v r (l0r t ) run number 











(M 


y- 1 ) 


(km s 4 ) 


in PSD 


CAK, FD case 
















0.6 


3.84 


X 


10~ 4 


1.9 x 


10 12 


10000 




0.6 


1.21 


X 


10~ 3 


1.3 x 


10~ n 


10000 




0.(5 


3.84 


X 


10~ 3 


8.7 x 


10~ n 


10000 




CAK, FD&Mmax, e case 
















0.4 


3.84 


X 


10~ 4 


2.0 x 


-i 17 

10 17 


700 




0.4 


6.03 


X 


10~ 4 


6.3 x 


10 ib 


1400 




0.4 


1.21 


X 


10 3 


6.3 x 


10 15 


3000 


15 


0.4 


3.84 


X 


10 3 


1.6 x 


10 13 


6000 


1(5 


0.4 


1.21 


X 


10~ 2 


2.4 x 


io- 12 


10000 




0.(5 


3.84 


X 


10" 4 


4.8 x 


1Q -14 


900 


2 


0.6 


6.03 


X 


10" 4 


5.5 x 


IO" 13 


2000 




0.6 


1.21 


X 


10~ 3 


4.7 x 


io~ 12 


3500 


3 


0.(5 


3.84 


X 


IO" 3 


4.0 x 


IO" 11 


5200^ 


4 


0.6 


1.21 


X 


lO" 2 


3.1 x 


10 -io 


7500 


5 


0.8 


3.84 


X 


10- 4 


1.9 x 


IO" 13 


500 




0.8 


6.03 


X 


10" 4 


8.8 x 


IO" 12 


1200 




0.8 


1.21 


X 


10~ 3 


6.2 x 


io- 11 


2500 


17 


0.8 


3.84 


X 


10~ 3 


6.3 x 


10 -io 


7500 


18 


0.8 


1.21 


X 


10- 2 


4.0 x 10~ 9 


14000 


19 



a) We re-estimated a typical velocity for this model. 



cal point. We can also express wind terminal velocity as a 
function of M(t c ): 



M max T the same as the fiducial model and the Eddington 
numbers as PSD's model 9. 



vl, oc ( * r ") = M(t c )T, 



(1 



(7) 



(see the Appendix). Thus luminous stars do not use all their 
available effective radiative force to drive winds. However 
for low luminosity stars, the importance of M max and the 
effective luminosity is essential as M ma x limits M(t c ). For 
example, equation (A22) shows that in a first order approxi- 
mation v^. oc M max T* (see equation 7). Physically, it means 
that if the force multiplier is saturated the radiation pres- 
sure accelerates the wind maximally which is proportional 
to M max Y,, in the single-scattering limit. Additionally, as 
we showed in Fig. 1 for the stellar case and PSD showed 
for the disc case, radiation pressure cannot drive a wind if 
MmaxF < 1//fd,c ^ 1.6, that is if the effective luminosity 
is less than the Eddington luminosity. 

To confirm further that disc winds behave in the same 
way as stellar winds we calculated a test model to compare 
with PSD's fiducial steady model for T t =T D = 1.21 x 1CT 3 
(their model 8) . For the test model all parameters are as for 
the fiducial model except for r_o = T* = 3.84 x 10 -3 , in- 
creased by a factor of 10/7T, and M ma x = 1380 decreased by 
the same factor. With these parameters our test model has 



The results of our test model confirm expectations: the 
test model has M as PSD's model 9 with the same F and 
v r (r ) as PSD's model 8^ with the same FMmax- The test 
model also settles into a steady state. Then the disk mass- 
loss rate, as the stellar mass-loss rate, depends on the Ed- 
dington number while velocity depends on the effective lu- 
minosity at the critical point. 



3.2 Dependence of wind properties on a and k 



In the stellar case for M(t c ) < M max , the analytic solutions 
provide scaling of the wind properties on all model parame- 
ters. For example, the a parameter determines a slope of the 
power- law dependence of M, on luminosity and on k (i.e., 
M, oc (fcr«) 1/Q ). 

For M(t c ) > M m ax in the CAK, FDk,M max case, a 
and k determine the slope of a linear dependence of M* on 



t Note a typo in PSD's table 2: for model 8, M D should be 2.1 X 
IO" 15 Mq yr" 1 instead of 1.2 X 10~ 15 M Q yr" 1 . 
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Figure 2. Model mass-loss rates as function of the Eddington number for k = 0.2 (upper panel) and k rescaled with a as described in 
Section 3.2 (lower panel). The thick lines represents results in the stellar CAK&zFD case for a = 0.4,0.6 and 0.8; the solid, dotted, and 
dashed lines, respectively. The thin lines represents results in the stellar CAK, FDlkM max case, the line types are as in the C'AK&zFD 
case. The open circles represent results for the disc CAK, FD&zM m a X ,e case and a = 0.6, while the crosses and asterisks represent results 
for a = 0.4 and a = 0.8, respectively. The alternative ordinate on the right hand side of the two panels is as in Fig. 1(a). 



luminosity: M* cx Fk 1 ^" . Additionally, a enters the expres- 
sion for M(t c ) and then determines the lowest value of F for 
which M(t c ) < M max and the power-law applies (see equa- 
tion 5). Moreover, through /fd,c the a parameter defines 
the lower limit of F for which a wind solution exists (see 
equation 6). 

Fig. 2 illustrates how the mass-loss rate depends on 
the Eddington number for various a with (a) k = 0.2 and 
(b) k rescaled with a (see below). Fig. 2 also compares 
our predictions for M* and Md- We do not consider the 
CAK, FD, &zM max case with a = 0.4 and 0.8 because our 
calculations for a — 0.6 are sufficient to show the basic be- 
havior in this case. 

In the disc case, we used PSD's results complemented 
by our results. PSD calculated a few disc wind models with 
a = 0.4, 0.6 and 0.8 for various disc and stellar luminosities 



in the CAK, FD&zM max ,e case. They found that there is no 
change in the power-law dependence of Md on luminosity 
between a — 0.6 and a — 0.8, in contrast to scaling found 
for the stellar wind solutions. 

The analysis in previous subsection showed that Md 
scales with the Eddington number for a = 0.6 almost as 
M» . However the analysis also showed that the CAK scaling 
does not apply in both the disc and stellar cases if T is 
very low and a finite M max is used. Using these conclusions 
about similarities between disc and stellar winds, a simple 
interpretation of the PSD results for a = 0.8 is that for 
these disc models M(t c ) < M max and the power-law does 
not apply. As we mentioned above, M(t c ) increases with a 
(equation A17) and subsequently the lowest F for which the 
power- law applies, increases to unity with a, (equation 1). 
Then for a fairly high a = 0.8, even the highest-considered 
r in Fig. 2 might still be in the region where the power-law 
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does not apply in the disc CAK, FD$zM max ,<, case. If so, no 
change in the power-law found by PSD and confirmed here 
can be seen as cancellation of the two opposite trends 1) a 
decrease of the slope caused by higher a and 2) an increase 
of the slope caused by a finite M max (see Fig. 1). 



that the mass-loss rate increases with decreasing a. This 
result can be understood if we consider the coupling be- 
tween the parameters of the force multiplier. As equation 
(9) shows, k increases with decreasing a causing an overall 
increase of M(t) and subsequently an increase of M. 



This interpretation is confirmed by the results for disc 
winds with a — 0.4. Fig. 2 shows that for a = 0.4, the slope 
is as predicted by the CAK&FD scaling and the power-law 
applies for a Eddington number range wider than for a = 0.6 
and a = 0.8. 

In their preliminary parameter survey, PSD calculated 
disc wind models for different a and T o while holding Mmax 
and k fixed. In their formalism, they did that by means of 
adjusting a value of r/ max according to the equation 



PSD's results suggest that M» increases with a. This 
result is the consequence of changing a but formally keeping 
k constant. We found a similar behaviour of M* (see Fig. 
2a) when we applied PSD's approach. However our previous 
analysis and numerical tests show that we can rescale M D 
if other k is chosen. For example, if we take PSD models for 
a = 0.4 and 0.8 and want to rescale them, keeping Mmax and 
rjmax as for a = 0.6 (i.e., M max = 4400 and r/ max = 10 7 ' 9 ) 
we need first calculate k n using equation (9) and then use 
the correction factor fk ■ In the example considered here the 
correction factor is: 



As we mentioned in Section 3.1 Gayley (1995) showed that 
M m ax can be interpreted as a pseudo-conserved quantity 
that is proportional to metallicity and is estimable a pri- 
ori. However k depends on metallicity, sound speed and the 
other force multiplier parameters (e.g., CAK, Gayley 1995). 
Assuming that Mmax is constant and using equation (4) we 
can express k as 



k — r/max Mmax- (9) 

1 — a 



Fig. 2b shows results from Fig. 2a, where the mass-loss rates 
for a — 0.4 and 0.8 are rescaled as described. The results 
for a — 0.6 arc unchanged. Fig. 2a clearly shows that after 
rescaling Mo and M* increase with decreasing a as found 
in previous stellar wind models. 

Contrary to the mass-loss rate, the terminal velocity 
does not depend on k (e.g., equations A13, A22, and A25). 
Our numerical tests confirm this for stellar and disk winds. 



Additionally adopting Gayley's working assumption that we 
truncate the line distribution at the single strongest line 
(i.e., M max = ^r-rjmax using our notation) we can rewrite 
equation (9) as 

^=T^S-fT M -- (l0) 

The last equation shows explicitly how the formal parameter 
fc depends on physical properties of wind. 

In this context we could ask how PSD's and our results 
for fixed k and M m ax can be rescaled to accommodate the 
coupling of k with other model parameters. 

Using the result for stellar winds: M oc k 1 ^", we can 
rescale the mass-loss rates for a new k parameter, k n if we 
multiply M by the correction factor: 

* = (t) • 

The detailed calculations of the force multiplier and 
stellar winds (e.g., CAK; Abbott 1982; Gayley 1995) showed 



3.3 Dependence of wind properties on sound 
speed 

Analytic steady state solutions of stellar winds depend very 
weakly on the sound speed (e.g., CAK; PPK). Numerical 
simulations of stellar winds showed that a higher c s short- 
ens the computational time needed for the initial conditions 
to settle into a steady state (e.g., OCR). However the final 
steady state solution does not depend on c s . Our numeri- 
cal tests confirm these results for stellar winds. We thus do 
not expect changes of our disc wind results with c s either. 
Nevertheless we calculated a few models for various c s . 

We have calculated models with c' s higher and lower 
than c' s — 4.6 x 10~ 3 by a factor of 3 and compare them with 
PSD's fiducial complex and steady models, PSD's model 2 
and 8 respectively. We found that Md and velocity of the 
disc winds do not depend on c' s . However, the time behaviour 
of the winds does depend. For example, the PSD fiducial 
steady state model recalculated with c' a = 1.53 x 10~ 3 is 
somewhat time dependent: density fluctuations originating 
in the wind base spread in the form of streams sweeping 
outward. We did not find infalls of dense clumps of gas on 
the disc. On the other hand, the model as the PSD fiducial 
complex model but with c' s = 1.38 x 10~ 2 settles into much 
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smoother and steadier state than PSD's model 2. Generally, 
we found that wind time dependence weakens with increas- 
ing c s because the gas pressure effects get stronger with 
increasing c' s . Subsequently higher gas pressure smooths the 
flow more effectively and in a larger region above the disc 
mid-plane as the size of the subcritical part of the flow in- 
creases with c' 



4 DISCUSSION 

Our comparison between disc and stellar winds has shown 
that despite differences in geometry, these two winds are 
similar in terms of M and «oo • In particular, we showed that 
M* is an upper limit on M D for the same model parameters, 
typically M„/M D <, 2. 

Previous CV wind simulations showed that the pre- 
dicted Md is lower than the observed Md by a factor of a 
few, at best (e.g., Drew 1997, PSD and references therein). 
If we use M» for the Eddington number and other model pa- 
rameters as for CV's accretion discs we still may not repro- 
duce observed disc mass-loss rates. Thus we might conclude 
that 1) present disc wind models still need some important 
improvement or 2) Md inferred from observations are over- 
estimated, or both. At the moment redetermination of Md 
from observations taking into account wind structure pre- 
dicted by PSD is probably better than moving directly to 
exploring different ways of enhancing theoretical Md (e.g., 
including magnetic field effects). Such an approach will pro- 
vide better constraints on our disk wind model. 

Our tests with different sound speeds showed that ther- 
modynamic effects might be important in determining time 
behavior of the disc winds. However we do not believe that 
a sound speed higher by a factor of 3 compared to the one 
used PSD for CVs is possible because this could correspond 
to the wind temperature higher by a factor of 9! Thus we do 
not expect great changes in models for CV disc winds after 
the inclusion of temperature calculations. 

Thermodynamics and changes in wind photoionization 
might be important in very luminous systems such as AGNs 
and X-ray binaries. In studies of the former there is the long 
standing problem of wind over-ionization by the central ra- 
diation, especially in the inner part of an outflow. Such over- 
ionization will manifest in a low M max , for instance. How- 
ever it is arguable how damaging for radiation driving that 
can be. Using fiducial values for AGN's: V = 0.4 and a = 0.6 
(e.g., Murray et al. 1995; Arav 1996) and equation (A19), 
M(t c ) = 3.75. Our analysis showed that line driven winds 
are strong and fast when M max > M(t c ). Thus for AGNs to 
produce strong and fast line-driven winds not too many ab- 
sorption lines are needed. Stevens & Kallman (1990) showed 
that M m ax < 4 for relatively high photoionization parame- 
ters (see also Arav, Li& Begelman 1994; Murray et al. 1995). 
This simple argument supports the idea that winds in AGNs 
are line-driven. 



To conclude we would like to mention that previous at- 
tempts to model disc winds using ID approximations, have 
found some vindication in our fully 2D time- dependent sim- 
ulations. This is not surprising, because of the energy argu- 
ment: in both the stellar and disc cases most of the radiation 
energy is released near the central object. Thus the fully 2D 
calculations are crucial to determine the density and veloc- 
ity as function of the position and time but less necessary 
to estimate the mass-loss rate and typical velocity for disc 
winds. 
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APPENDIX A: SUMMARY OF ANALYTIC RESULTS FOR THE CAK MODEL WITH THE FINITE 
DISC CORRECTION AND A FINITE MAXIMUM FORCE MULTIPLIER. 

The radial radiation force (per unit mass) due to spectral lines at every location in the flow can be expressed as 

F rad > l = J (£«f^M(t), (Al) 

where the term in brackets is the radiation force due to electron scattering from a stellar surface element and dT* is the 
frequency integrated flux from that element. The integration is over the visible radiant surface. The quantity M, called the 
force multiplier, represents the increase in radiation force over the pure electron scattering case when lines are included. CAK 
found a fit to their numerical results for the force multiplier, such that 

M(t) = kt~ a , (A2) 
where k and a are parameters of the fit. In the Sobolev approximation, M(t) is a function of the optical depth parameter 
dvi - 1 



t = a e pv th 



dl 



(A3) 



where p is the density, v t h is the thermal velocity, and -4 s - is the velocity gradient along the direction linking a point in the 
wind to the radiant surface element. For a spherically-expanding wind, 

dVl dv r 2 V r 2 X , . .n 

-dT = ^ + v {1 -^' (A4) 

where v r is the radial wind velocity and p is the direction cosine. 

Making an assumption ^ 3> ^- , certainly valid in a lower part of the wind acceleration zone, ^ is 

dvi dv r 2 , a ,n 

If = ~~d~F^ ' (A5) 

Then using equations (A3) and (A4), the equation for F ' is 

p ra d ,l = (^fl^ f FD M C AK(t), (A6) 



where is the frequency integrated flux from the star, McAi<(t) = k [a e pvth \ ^jjr \ j , is the force multiplier in the form 
introduced by CAK, and /fd is the so-called finite disc correction factor which allows a generalization of the CAK formalism 
beyond the point-star case. The assumption made above simplifies the finite disc factor to a function of radius only: 



l _(l_ r ;/ r 2 ^ 
(1 + a)r»/r 



fMr) = , J'2,'2 ■ (AT) 



Using the Eddington number, T* = 4 ^*q^ , where all symbols have their conventional meaning, we can rewrite equation 
(A6) as 

F rad ' L = 9MI±f FD {r)kt- a . (A8) 
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The equation of motion for a spherical, steady state wind is 

v _ S±) — + GM *( l - T *) + M _ M + i? rad - ; = o (A9) 
v J dr r 2 r dr 

where c s is the gas sound speed. We dropped the 'r' subscript for simplicity. The principles of solving this equation of motion 
are presented in CAK (see also PPK). Here we generally follow them to calculate the velocity law and the stellar mass-loss 
rate in the zero sound speed case. Equation (A9) reduces to 

^ + GM.(i-r.) + ^ l = 

dr r 2 

Using the continuity equation and equation (A8), equation (A10) becomes: 



, , GM. l-T, , GMr, . . / 4tt 

r \a e v th M, 



mjl + ^_ 1 K_ H + ^±± fFD(r)k ( ) („„')« = , (All) 



where w' = f . 

Assuming /fdM = 1, we have the case considered by CAK. Here we repeat their well known solutions for the velocity 

law 

v(r) = Uoo (l-r*/r) /3 . (A12) 
where /3 = 0.5, 



2 a f 2GM,(1 -r») 



t4 = i -^-^ (A13) 

1 — a \ r* J 

and the mass-loss rate 

m. = 47rGM * — (i - r,r (1 " a)/a ((i - ^fer,) 1 /" . (Ai4) 

o^th 1 — a 

When the star is treated as a point-source, the optical depth parameter does not vary with radius and is 

t = (A15) 
and subsequently 

On the other hand using the critical point conditions (see CAK, PPK and Gayley 1995) the force multiplier is 

MW=M(ic) = (r^k' (A17) 
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where the 'c' subscript indicates that a quantity is evaluated at the wind critical point. 

Comparing the right hand sides of equations (A16) and (A17) and using equation (A13), M* can be expressed as 

^ = 4^M. ( i-r.)_w k y 



CT e Vth 1 — OL \M(t c 

Equation (A18) demonstrates the importance of the force multiplier at the critical point in determining the mass-loss rate. 

For the finite disc case, the mass-loss rate is also given by equation (A18) but the force multiplier at the critical points 
is now 

M(Q = 1 ~ T * . (A19) 

(1 — a)l ,jFD,c 

For the velocity law we assume that it can be approximated by equation (A12) and find that 

vl = _ r . /re) -* + i. (A20 ) 

We confirm PPK results that this approximation works fairly well for (5 — 0.8 and r c — 1.05r*. Note that in this case the 
force multiplier increases with radius. 

Now we move to a problem of solving the equation of motion when M max < M(t c ). In this case we assume that the force 
multiplier is constant, i.e., M(t) = M max for all locations in a wind. Such a simple approach allows the equation of motion to 
be written as 

, , GM,(l-r,) GMr, , n , ao1 , 

w + ^ '- + —^—f F D(r)M max = 0. (A21) 

For /fd = 1, we find that the velocity law is again as in the CAK case with 
^ = 2GM.(r.(l + M ma3! )-l) | 

while the mass-loss rate is 
^ _^GM,(Y,{1 + M max )-1) ( k \ lf 



iw-) • (A23) 



For the finite disc case, the velocity law is: 

v ( r ) = ((iTST^ 1 +«)(!- r*) r * + r,M max (r 2 (l - H2Fl{-0.5, -a; 0.5; $)) 

V 2 xi/2 (A24) 

-r 2 H2Fl(0.5, -a; 1.5; £)) + rr.(r.(l - M max (l - hi - h2)) - a(l - r,) - 1))J 

where H2F1 is hypergeometric function, hi = H2F1(— 0.5, —a; 0.5; 1) and h2 = ff2Fl(0.5, —a; 1.5; 1). The terminal velocity 
is 
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2GM, 



1/2 



v.. 



(1 + a)r* 



- M max (l + a/r, -hi- h2)) - a{\ - r„) - 1) 



(A25) 



And finally the mass-loss rate is 



M, = 



a-kGM* (r. (i + f FD , c M max ) - 1) 



0~eV t h 




(A26) 
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